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We report several results concerning W(A, q) = exp(So/kB), the exponent 
of the ground state entropy of the Potts antiferromagnet on a lattice A. First, 
we improve our previous rigorous lower bound on W(hc, q) for the honeycomb 
(he) lattice and find that it is extremely accurate; it agrees to the first eleven 
terms with the large-g series for W(hc, q). Second, we investigate the het- 
eropolygonal Archimedean 4 • 8 2 lattice, derive a rigorous lower bound, on 
W(A ■ 8 2 ,q), and calculate the large-g series for this function to 0(y 12 ) where 
y = l/(q — 1). Remarkably, these agree exactly to all thirteen terms calcu- 
lated. We also report Monte Carlo measurements, and find that these are very 
close to our lower bound and series. Third, we study the effect of non-nearest- 
neighbor couplings, focusing on the square lattice with next-nearest-neighbor 
bonds. 
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I. INTRODUCTION 



Nonzero ground state disorder and associated entropy, So 7^ 0, is an important subject 
in statistical mechanics; a physical realization is provided by ice, for which So = 0.82 ± 0.05 
cal/(K-mole), i.e., S /k B = 0.41 ± 0.03 0. A particularly simple model exhibiting ground 
state entropy without the complication of frustration is the g-state Potts antiferromagnet 
(AF) on a lattice A for sufficiently large q p|- 0. This subject also has a deep connection 
with graph theory in mathematics @|- since the zero-temperature partition function 
of the g-state Potts AF on a lattice A satisfies Z(A,q,T = 0)paf = -P(A, q) and hence 
So/ks = lnW(A,q), where P(G,q) is the chromatic polynomial, expressing the number of 
ways of coloring the vertices of a graph G with q colors such that no two adjacent vertices 
have the same color, and 

W(A,q) = limP(A n ,g) 1 / n (1.1) 

where A n denotes an n- vertex lattice of type A (with appropriate boundary conditions). 
Given the above connection, it is convenient to express our bounds on the ground state 
entropy in terms of its exponent, W(A,q). Recently, we studied the ground state entropy 
in antiferromagnetic Potts models on various lattices, including both Monte Carlo measure- 
ments and rigorous upper and lower bounds J7]- JTD|. Here we continue this study. The 



reader is referred to Refs. M- |1 for further background and references. 



II. IMPROVED LOWER BOUND ON W(hc, q) FOR THE HONEYCOMB LATTICE 

In Ref. 0, we derived lower and upper bounds on W(hc, q) for the honeycomb (he) 
lattice using a method first applied by Biggs to obtain such bounds for the square lattice 
| i~lfl . We showed that the upper and lower bounds rapidly approached each other for large 
q and hence became very restrictive even for moderately large q. In particular, we obtained 
the lower bound, W(hc, q) > W(hc, q)e> for q > 3, where 

(a - I) 3 / 2 

W(hc, q ) e = {q ql / 2 (2.1) 

and found that this was very close to the actual value of W(hc, q), as determined from Monte 
Carlo measurements 10-0, [fT2| . (We also obtained similar bounds for the triangular lattice.) 
Accordingly, in this paper we shall focus on deriving rigorous lower bounds for several lattices 
A, which again turn out to be very accurate not just as bounds but as approximations for the 
W(A, q) functions. In passing, we mention that, using the same methods as in Ref. 0, we 
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could also derive rigorous upper bounds for additional lattices; however, we concentrate here 
on the lower bounds because of their very high accuracy as approximations for the W(A, q) 
functions. We begin by sharpening our lower bound for the honeycomb lattice. We picture 
this lattice as a brick lattice with the long axis of the bricks horizontal, as in Fig. 1(a) of 
Ref. ||. We consider a sequence of these lattices with m (n) sites in the horizontal (vertical) 
direction, denoted (hc) mxn . In the thermodynamic limit the boundary conditions do not 
affect the bound, so for technical convenience, we take free (periodic) boundary conditions 
along the horizontal (vertical) directions. Further, to maintain the bipartite nature of the 
lattice and avoid frustration, we take n even. Our bounds in Ref. Q were obtained by 
taking the vectors of color configurations to refer to the vertical lines; here we take them 
to refer to horizontal lines L. The number of allowed g-colorings of this line of vertices is 
M = P(L m , q) = q(q — l) m_1 . One can then associate with two adjacent horizontal lines of 
sites L, V , forming one layer of bricks, an M X M dimensional symmetric coloring matrix T 
with entries T^li = 1 or if the g-colorings of these lines are or are not compatible. Here, 
by compatible g-colorings we mean colorings that satisfy the constraint that no two adjacent 
vertices have the same color. Then (for fixed m,n) P((hc) mxn ,q) = Tr(T n ). Since T is a 
nonnegative matrix, one can apply the Perron-Frobenius theorem [13[] to conclude that T 
has a real positive maximal eigenvalue A maXjn (g). Hence, for fixed m, 

JimTr(n 1/(mn) ^Al- (2.2) 

so that 

W(A,q)=lim o \llZ (2.3) 

Denote the column sum Kj(T) = YhLiT^ (equal to the row sum Pj(T) = YnLiTji since 
T T = T) and S(T) = E^=i% note that S(T)/Af is the average row (column) sum. 
Combining the bounds for a general nonnegative M x jV matrix A, |]T3 



min{7 i (A)} < \ max {A) < max{7 j (A)} , for ^ = Kj or pj (2.4) 
with the {y = 1 case of the ) more restrictive lower bound applicable to a symmetric non- 



negative matrix 14 



M 
we have 

SiT) 



S(T U )~\ 1 / u 

' < \ max , for v = 1,2,... (2.5) 



< ^max{T) < max{« i (T)} (2.6) 



2 



The sum S(T) is given by the chromatic polynomial for the brick layer of length m, hence 
containing 2m vertices, which we calculate to be 

P((Ch) 6!2m , q) = q{q - l)D & {q)^-^ (2.7) 

where the notation JK| (Ch)k,2m means a 2m-vertex chain of fc-sided polygons, with each 
adjacent pair of /c-gons intersecting along one common edge, and 

s=0 V s I 



Taking the m — > oo limit and applying the v = 1 case of (2J3), we thus obtain the lower 
bound W(hc, q) > W(hc, q)i for q > 3 fl2f , where 

W(hc , q)t = EMU = ^-5^ + l<V-10, + 5)V> 
q — 1 q — 1 



For q > 4, the difference between the lower bounds ( |2.1| ) and (|2.9| ) decreases rapidly toward 
zero. The corresponding reduced W functions W r (hc,q) = q^Wihc^q) have large-g Taylor 
series which coincide up to order 0(g~ 4 ): 

3 3 1 3 131 

q^W(hc, q) t = l- 6 -q-' + |<T 2 + -<T 3 + ^g" 4 + ^g~ 5 + 0(g~ 6 ) (2.10) 

and 

q- X (W{^,q) t -W{h^q) e ) =\q- 5 + 0(q- e ) (2.11) 

Since the quantity y = l/(g — 1) is a natural variable for a large-g expansion, it is also useful 
to compare an expansion of our new lower bound with a small-?/ series. For a lattice A, this 



series has the form |Tj|- [I7| 

W(A, q) = q(l - q- r f 2 W(A, y) (2.12) 
where ( is the lattice coordination number, and 

OO 1 

W{A,y) = l + J2 w ny n , V= 7 (2-13) 

Analogously, for the expansion of our lower bound, we define the reduced lower bound 
function W(A, y)i via 

W(A, q) t = q(l - q-^W^A, y) t (2.14) 
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From eq. (|2.9|) , we find the very simple result 

W(hc,y) e = (l + y 5 ) 1 / 2 
Expanding this in a Taylor series in y, we have 

W(hc, V ) t = l + \y 5 - ^y 10 + ^y 15 + 0(y 20 ) 



(2.15) 



(2.16) 



Although as a lower bound, this need not, a priori, agree with terms in the small-?/ Taylor 
series expansion of the actual function W(hc,y), we find that, remarkably, it does agree all 
the way up to 0(y 10 ), i.e., for the first eleven terms. The small-y Taylor series expansion of 
W(hc, y), calculated to order 0(y 5 ) in Ref. |16j and to 0(y 18 ) in Ref. |L7j 



is: 



W(hc,y) e = l + ~y 5 -~y 10 + y n 



0(y 



V2\ 



(2.17) 



To show how close our lower bound is to the actual function W(hc, q), we give a compar- 
ison in Table 1 for q = 3 through q = 10. We determined the values of W(hc, q) by Monte 
Carlo measurements, as discussed in detail in our previous Refs. and ||. We also in- 
clude a comparison with our slightly less stringent lower bound fl2.1|) . Evidently, both lower 
bounds are quite close to the actual measurements, even for q as low as 3. In particular, the 
sharper lower bound (|2.9| ) provides an extremely good approximation to the actual function 
W(hc, q) for q > 3. 

TABLE I. Comparison of lower bounds for W(hc, q) and ratios of bounds to Monte Carlo mea- 
surements for 3 < q < 10. The estimated uncertainty from the MC measurements for the entries 
in the W(hc,q)i/W(hc,q)MC an d W (he, q)i> /W (he, q) mc entries is 3 x 10~ 4 . 



Q 


W(hc, q) M c 


W(hc,q) e 


W(hc,q) e 
W(hc,q) M C 


W(hc, q)i> 


W{hc,q) e , 
W(hc,q) MC 


3 


1.6600(5) 


1.6583124 


0.99898 


1.6329932 


0.98373 


4 


2.6038(7) 


2.6034166 


0.99985 


2.5980762 


0.99780 


5 


3.5796(10) 


3.5794553 


0.99996 


3.5777088 


0.99947 


6 


4.5654(15) 


4.5650849 


0.99993 


4.56435465 


0.99977 


7 


5.5556(17) 


5.5552778 


0.99994 


5.5549206 


0.99988 


8 


6.5479(20) 


6.5480952 


1.00003 


6.5479004 


1.00000 


9 


7.5424(22) 


7.5425874 


1.00002 


7.5424723 


1.00001 



10 8.5386(25) 8.5382220 0.99996 8.5381497 0.99995 
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III. 4 8 2 LATTICE 



Another interesting question in statistical mechanics and graph theory is the Potts AF 
ground state entropy and the related asymptotic behavior of the chromatic polynomial on 
heteropolygonal lattices, i.e. Archimedean lattices [18] involving tiling of the plane by more 
than one regular polygon such that all vertices are equivalent. In standard mathematical 
notation ||T8"1 , |T9| , an Archimedean lattice is defined by the set p\ ■ P2 • ■ ■ p n of polygons that 
one traverses in a circuit around a vertex (site) of the lattice. For a homopolygonal lattice, 
these are all the same (e.g., in this notation, the square, honeycomb, and triangular lattices 
are 4 4 , 6 3 , and 3 6 ). We consider here the 4 • 8 2 lattice. 



A. Lower Bound on W{A 



It is convenient to represent the 4 • 8 2 lattice as a generalized brick lattice, contain- 
ing two types of bricks at each horizontal level: square bricks and rectangular ones with 
(length, height) = (3, 1). In the next layer up, a square brick is located in the center above 
a (3,1) brick, and a (3,1) brick is located with its center above a square brick, and so 
forth for successive levels. Using the same general method as before, we construct a color- 
ing matrix defined between adjacent horizontal rows of bricks and obtain the lower bound 
W(4 • 8 2 , q) > W(A ■ 8 2 , q) t for q > 3 [0, where 



W(4-8 2 ,q) i 



D4(q)D 8 (q) 



1/4 



{q 2 -3q + 3)(g 6 - 7q 5 + 21q 4 - 35g 3 + 35g 2 - 21q + 7) 



1/4 



q-1 



Dividing by the asymptotic behavior as before |§, this has the large-g expansion 
q- 1 W{A-S 2 ,q) l 

z, o 

The corresponding reduced function W(4 ■ 8 2 ,y) is 



W(4-8 2 ,y) e 



[l + y 3 + y 7 + y 



l + y) 1/2 (l-y + y 2 )(l-y + y 2 



10^1/4 



1/4 



y 3 + y 4 



y 5 + y 6 ] 



1/4 



(3.1) 



(3.2) 



(3.3) 
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B. Large-g Series Expansion 



Using standard methods [ 15l , we have calculated the large q expansion, i.e., the small-?/ 
expansion, of the reduced function W(A ■ 8 2 , y). We find 

W(4 ■ 8 2 , y) = 1 + \f - |y 6 + \y 7 + ^ + ly 10 - ^ + 0(y 13 ) (3.4) 

Next, in order to assess the accuracy for large q of our lower bound, we have expanded 
the corresponding reduced lower bound function W(4 ■ 8 2 ,y)e, eq. ( |3.3|) in a series in y. 
Remarkably, the small-y series expansion of our lower bound agrees to 0(y 12 ), i.e. to all 
thirteen terms that we have calculated, with the small y expansion of the true W(4-8 2 ,y) 
function in (|3.4|) ! This shows that for even moderately large q, our lower bound IT(4- 8 2 , q)t 
actually provides an extremely accurate approximation to the true function W(4 ■ 8 2 , q). Of 
course, it is straightforward to calculate the series ( |3.4j) to higher order in y, and at 0(y 13 ) 
or above the series for W(4 ■ 8 2 , y) and W(4 ■ 8 2 , y)e could differ; nevertheless, the agreement 
to 0(y 12 ) is already quite striking. 



C. Monte Carlo Measurements 

To get further information on W(4 ■ 8 2 , q), we have carried out Monte Carlo (MC) mea- 
surements of the Potts AF ground state entropy So and hence obtained W = exp(S'o//cs), 
for 3 < q < 10. Our methods are similar to those used in our previous works 0- 0. We 
present our results in Table 2, including comparisons with our rigorous lower bound and 
our series calculation to 0(y 12 ) (to which order they coincide). Evidently, the agreement is 
extremely good, even for q as low as 3. 

TABLE II. Exponential of ground state entropy, from Monte Carlo measurements, and compar- 
ison with lower bound, for the Potts antiferromagnet with 4 < q < 10 on the 4 • 8 2 lattice. The 
estimated uncertainty from the MC measurements for the entries in the W(4- 8 2 , q)i/W(4:-8 2 , q)MC 
column is 3 x 10~ 4 . 

* W(4-8\q)MC gggafc 

3 1.68575(60) 0.99959 

4 2.62226(75) 0.99994 
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5 3.5918(10) 0.99995 

6 4.5737(14) 0.99996 

7 5.5618(16) 0.99992 

8 6.5530(20) 0.99995 

9 7.5466(23) 0.99995 

10 8.5413(26) 0.99998 
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IV. SQUARE LATTICE WITH DIAGONAL NEXT-NEAREST-NEIGHBOR 

COUPLINGS 

An interesting question concerns the effect of non-nearest-neighbor interactions on the 
ground-state entropy of the Potts antiferromagnet and the corresponding VT(A, q) function. 
We explore this question here by considering the Potts AF on a square lattice with next- 
nearest-neighbor couplings along both diagonals, defined by the Hamiltonian 

W = -E^ rf (4-1) 

n, n' 

where J n y = J < if the sites n and n' are nearest neighbors along a row or column, 
and J n>n i = Jd < if n and n' are diagonally opposite sites on a given square of the lattice 
(with each such pair counted only once in the sum (4.1)). We further define K = f3J and 



Kd = (3Jd and consider the T — ► limit, so that the color on each site n is required to be 
different from the colors of each of its nearest neighbors on the rows and columns of the 
lattice, and also each of its next-nearest neighbors along the diagonals of the lattice. Then 
Z((sq)d,q,T = 0)paf = P((sq)d, q), where (sq)d denotes the nonplanar graph formed from 



a square lattice by adjoining bonds connecting diagonally opposite sites on each square [pO 
Note that the (sq)d lattice has coordination number £ = 8 and girth 7 = 3, where "girth" 
7 of a graph G denotes the length of the shortest circuit on G. We find that the chromatic 
number of this lattice (for free boundary conditions or periodic boundary conditions that 
preserve the bipartite structure of the underlying square lattice before inclusion of diagonal 
bonds) 

X((sq)d) = 4 (4.2) 

(x(G) is defined in footnote [0.) These are to be compared with ( = 7 = 4 and x( S( l) — 2 
for the square lattice. From our earlier discussion, one thus has W((sq)d,q) < W(sq, q). 
Indeed, we find that there is an especially strong difference with respect to the square lattice 
for the case q = x(( S( l)d) = 4. Here, we calculate that for a (sq)d lattice A LlxL2 of length L\ 
in the x direction and L 2 in the y direction, 

P(A LlxL2 ,4)~2 Ll+i2 (4.3) 

in the thermodynamic limit. Although this diverges, it does not diverge rapidly enough to 
yield a finite entropy. With the ordering of limits 

W(A,q s ) = lim lim P(G,q) 1/n (4.4) 

n— >oo q^q a 
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as discussed in Ref. |§, we thus obtain the exact result 

W((sq) d ,q = 4) = l i.e., S ((sq) dj q = 4) = (4.5) 

This may be compared with the value for the usual square lattice, for which W(sq, 4) = 
2.3370(7) ||, corresponding to the nonzero ground state entropy So(sq, 4) = 0.8489(3). 
(Indeed, even without any detailed calculation, one knows that So(sq,4) is nonzero from the 
elementary rigorous lower bound on a bipartite lattice, So(Am p ., q) > (1/2) ln(g — 1).) 

TABLE III. Exponential of ground state entropy, from Monte Carlo measurements, and com- 
parison with lower bound, for the Potts antiferromagnet with 5 < q < 10 on the square lattice with 
diagonal next-nearest-neighbor couplings. The estimated uncertainty from the MC measurements 
in the entries for W(sq d , q)e/W(sq d , q)MC is 1 x 10 -3 . 



Q 


W(sq d ,q) M c 


W(sq d ,q) t 
W{sq d ,q)MC 


W(sq,q) M c 


5 


1.5781(16) 


0.95051 


3.2510(10) 


6 


2.4460(24) 


0.98119 


4.2003(12) 


7 


3.3660(33) 


0.99029 


5.1669(15) 


8 


4.3093(43) 


0.99453 


6.1431(20) 


9 


5.2680(53) 


0.99658 


7.1254(22) 


10 


6.2363(62) 


0.99774 


8.1122(25) 
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Using the same type of coloring matrix method as above, we obtain the rigorous lower 
bound W((sq) d ,q) > W((sq) d ,q)t, where 

W((sq) d ,q) e = iq ~ 2)iq ~ 3) for q>A (4.6) 
q - 1 

For large q, we have the Taylor series expansion 

q^Wisqa, q)t = l- Aq~ l + 2{q~ 2 + q' 3 + q~ 4 + 0(<T 5 )) (4.7) 
From eq. (|2.14j) we have 

W(sq d , y) e = (l-y)(l-2y)(l+y) 3 (4.8) 

The MC measurements are given in Table 1. For comparison, we include also our MC 
measurements of W(sq, q) ||. The lower bound ( |4.6|) steadily approaches the measured 
values W((sq) d , q) as q increases and is quite close to them for q greater than about 6. The 
inequality W((sq) d , q) < W(sq,q) is also evident. 

This research was supported in part by the NSF grant PHY-93-09888. 
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